Density of States Driven Anisotropics Induced by Momentum Decoupling in High-T c 
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Momentum decoupling arises when small momentum transfer processes dominate the electron- 
phonon scattering and implies that anisotropies in superconductivity are driven by the anisotropies 
of the density of states. Considering an isotropic s-wave interaction in the momentum decoupling 
regime we give a natural simultaneous explanation to various aspects of ARPES and tunnel exper- 
iments on BiiSr-zCaCuiOzi including the correlation of gap magnitude and visibility of the dip 
above the gap, the enhancement of anisotropy with temperature, the presence of gap minima away 
from the F — X direction and a gap maximum in the T — X direction, the similarity of tunnel and 
ARPES spectra in the Y — M direction and the asymmetry in the SIN tunnel spectra where the 
dip structure is present only at negative sample bias. 



The observation by Angle Resolved Photoemission 
Spectroscopy (ARPES) of the superconducting gap in 
BiiSr^C aCuiO% ||,[l],|],[I|||] failed in answering the con- 
troversial question of the symmetry of the order param- 
eter. Some experimental results support the d-wave hy- 
pothesis , while others point to a rather complex mixed 
state in which the gap has nodes away from the T — X 
direction Q or even has no nodes at all HQ ■ These con- 
tradictions, together with the absence of gap observation 
in Y Ba^CuzO-; and the weakness of the gap values re- 
ported on Bi 2 Sr2CaCu20s [[§, illustrate the difficulty in 
measuring the gap with the ARPES technique. 

However there are some qualitative points on which all 
ARPES experiments are in agreement. The first point is 
that the higher values of the gap (which are the easier to 
measure) are reached in the direction where the density 
of electronic states (DOS) on the Fermi level is maximal, 
and in general the gap value is smaller when the DOS 
is smaller. The second important point is that the dip 
structure above the gap is more visible in the direction 
where the gap and the DOS are maximal, and that the 
visibility of the dip structure follows the variations of the 
gap and DOS. A third important remark is that from the 
ARPES experiments in the optimal direction for the gap, 
one obtains a spectral function very similar to the tunnel 
spectra ||. In the case of tunneling we see the average 
of the spectral function over the Fermi surface, and the 
similarity of ARPES and tunnel data in so anisotropic 
materials is quite surprising. In addition, some very im- 
portant experimental trends have been reported in recent 
ARPES and tunnel experiments. In Ref. |dj it is shown 
that the anisotropy is strongly enhanced when we move 
from the T — regime to the T — > T c regime. In Ref. 
H it has been shown that the gap has minima about 
10° away from the T — X direction and a maximum in 
the r — X direction. On the other hand, detailed vac- 
uum tunneling spectroscopy measurements report an 
asymmetric density of states in Bi^Sr^C 'aC^Og+s, with 
the dip appearing only at negative sample bias. 

In the following, we will give a simultaneous expla- 



nation to all the previously cited experimental points. 
We will see in particular that all the previous points 
characterize a superconductor in which the isotropic s- 
wave electron-phonon (or other boson) coupling is domi- 
nated by forward scattering processes. Dominance of for- 
ward scattering processes can be due to the vicinity of a 
phase separation instability ||, that appears reasonable 
in models describing strongly correlated electronic sys- 
tems On the other hand, there is evidence from Ra- 
man scattering [Io| that small momentum transfer (Ra- 
man active) phonons are strongly affected by supercon- 
ductivity in cuprates. Notice that a somewhat analogous 
situation has already been considered by Abrikosov in 
order to understand the change of sign of the order pa- 
rameter reported in some experiments on YBCL2CU3O7 
pr[ , and that dominance of forward scattering is an opti- 
mal condition for a positive contribution of non-adiabatic 
effects to superconductivity ]E^ , |r3[ . 

When forward scattering is dominant, there is "Mo- 
mentum Decoupling" (MD) in the superconducting be- 
havior, implying a different coupling in different regions 
of the Fermi surface. Indeed one has to keep in mind 
that the coupling is an intensive quantity and not an 
extensive one. In the case of MD the coupling at each 
region of the Fermi surface is proportional to the local 
DOS, and therefore the anisotropies in the superconduct- 
ing state are induced by the anisotropies of the density of 
states in the normal state. 

Let us now illustrate briefly how the MD appears for 
small momentum transfer processes and why it is the only 
situation which leads to DOS dependent anisotropies. 
The anisotropic Eliashberg equation in the off-diagonal 
sector, for an Einstein spectrum can be written with 
usual notations as follows. 
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The k dependence is contained in the coupling \g(k— p)\ 2 ■ 
In conventional s-wave superconductors it is assumed 
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that the interaction \g(k — p)\ 2 is almost constant on the 
Fermi surface and it leads to an isotropic On the other 
hand, if one supposes that \g(k — p)\ 2 has a relevant mo- 
mentum dependence in the vicinity of the Fermi surface 
(as it is the case in the d-wave scenario where this func- 
tion reflects electron-spin fluctuation coupling) then we 
obtain a fc-dependent coupling and from equation (1) one 
can obtain an anisotropic gap. However the anisotropy 
of the superconducting parameters is mainly imposed 
by the anisotropy of the interaction and not from the 
anisotropy of the density of states. In order to obtain sig- 
nificant DOS induced anisotropics one has to consider an 
isotropic s-wave interaction dominated by forward scat- 
tering processes. This can be illustrated by taking for 
example an interaction which is sharply peaked at zero 
momentum \g(k — p)\ 2 ~ g 2 S(k — p). Then from equation 
(1), it is easy to see that there is momentum decoupling. 
We obtain a momentum independent Eliashberg equa- 
tion which provides the gap function A(k,iuj n ) for each 
momentum k on the Fermi surface. This last equation 
is analogous to the isotropic Eliashberg equation with a 
coupling strength proportional to the density of states at 
the given point of the Fermi surface N(Ep, k). 

Of course a ^-function peak at q — is a rather un- 
realistic coupling function. However MD occurs even for 
finite q provided q is small compared to the characteris- 
tic momentum of the DOS variations over the Brillouin 
zone. To illustrate this point we performed numerical 
simulations on a simple two-dimensional BCS model [fli) . 
In that case the gap is given by 

(2) 

We consider an isotropic s-wave electron-phonon cou- 
pling having at small momenta a Lorentzian behavior 
as a function of the norm of the exchanged momentum 
V(q) = —V(l + \q] 2 /q 2 ) In this spectrum the electron- 
phonon scattering is dominated by the processes which 
transfer a momentum smaller than q c . 

For clarity, we will consider here |l4f| the simple near- 
est neighbor tight binding dispersion at half-filling ^ = 
—2t[cos(k x ) +cos(k y )] (the lattice spacing is taken equal 
to unity). The Fermi surface is a square defined by 
k x = k y ± 7r and k x = —k y ±ir with saddle points at 
(0,±7r) and (±7r,0). The minimum of the density of 
states is obtained at the points (±7r/2, ±7r/2) and there- 
fore the characteristic length of the DOS variations over 
the Brillouin zone is ir/\/2. We expect therefore that for 
q c > the gap might be isotropic while for q c suf- 

ficiently smaller than MD should manifest leading 

to DOS induced anisotropics. In fact in figure 1 we show 
the ratio of the gap at (0, it) over the gap at the points 
where the DOS is minimal (7r/2, 7t/2) as a function of q c . 
We can see that for q c < this ratio begins to be 



appreciably different from unity indicating the onset of a 
DOS induced anisotropy because of MD. 

We are now going to see how MD can explain all the 
features of the ARPES and tunnel experiments men- 
tioned in the introduction. We take first the tempera- 
ture dependence of the anisotropy reported in Q . If MD 
was perfect in BiiSriCaCuiO^+i, the temperature at 
which the gap disappears in the Y — X direction should 
be smaller to that in the Y — M direction. In fact since 
the DOS is smaller in this direction, the coupling and T c 
should also be smaller. Therefore, in the case of MD the 
anisotropy is enhanced close to T c , and if MD would be 
perfect it should even diverge. 

We will now show how within our BCS model the re- 
sults of § can be qualitatively reproduced when finite 
(small) momenta are transfered. In figure 2 we report 
the temperature dependence of A(0, ir) and A(ir, tt) for 
different values of q Cl and in figure 3 we give the corre- 
sponding temperature dependence of the anisotropy ratio 
R = A(0, 7r)/A(7r, tt). The critical temperatures are ob- 
tained by solving numerically the hermitian eigenvalue 
problem of the linearized equations near T c |l4| . 

When q c = 7r/4 the DOS induced anisotropy, because 
of partial MD, is already significant (R sa 1.7) but the 
anisotropy is almost temperature independent (Fig. 3). 
For smaller values of q c there is a continuous deforma- 
tion of the T dependence of A(7r, tt) |l4| from the large 
q c regime to the q c — > regime where, as expected in 
perfect MD, A(tt.tt) should have a BCS behavior going 
to zero at a temperature of the order T c /2, and therefore 
the anisotropy ratio should diverge close to T c . The re- 
sults of Ref. Q point to a strong MD regime. But one 
should bear in mind that, if the small gap is smaller than 
the temperature at which it is measured, it becomes ex- 
perimentally inaccessible flBl . Taking into account these 
damping effects neglected in our BCS model, the results 
of Ref. [|4| can be qualitatively understood even with q c of 
the order 7r/10 |14f . The enhancement of anisotropy with 
temperature is an evidence of MD, and cannot be under- 
stood in the case of anisotropic interactions like those 
considered in the d-wave scenario. 

It is difficult to conclude on the value of the coupling 
from the gap measurements made by ARPES, because of 
their large uncertainties. However there is a qualitative 
feature common to all experiments that certifies that in- 
deed moving from the T— M to the Y—X direction on the 
Fermi surface we go from a strong coupling regime to a 
weak coupling regime. In fact in the Y — M direction one 
can observe a dip structure above the gap that is a strong 
coupling effect independent on the spectral structure of 
the boson Such a dip appears when for sufficiently 
strong couplings the gap A is comparable to the boson 
energies that mediate superconductivity (2A/T C > 5.5 ). 
The stronger is the coupling, the sharper and deeper is 
the dip Q. It is a common trend of all experiments that, 
moving from r — X to Y — M , the visibility of the dip 
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follows the enhancement of the gap and DOS, indicating 
the presence of different couplings at different regions of 
the Fermi surface in agreement with the MD picture. 

In Ref . , it has been reported that the gap has min- 
ima around 30° — 35° and 60° — 65° (angles measured 
from X, Y — M) and a smaller maximum at 45°. The au- 
thors claim that there must be a difference between the 
results in the TX and TY quadrants (which is not ap- 
parent from their data) and that the minima might cor- 
respond to nodes. To our analysis instead, these results 
will be shown to be an evidence that gap anisotropies are 
indeed driven by DOS anisotropies, indicating that the 
interaction is isotropic s-wave in the MD regime. 

We show in figure 4a, the angular dependence of the 
density of states (measuring angles exactly as in Ref. ) 
for the simplest next nearest neighbor tight binding dis- 
persion that accounts qualitatively for the CuO bands 
seen by ARPES @. Fixing the distance fr om the bot- 
tom of the band at ss 350meV as in the experiment || , 
we consider three different characteristic situations de- 
pending on the distance of the van-Hove singularity from 
the Fermi level SE (SE = lOmeV , AQmeV and 90meV). 
We can see that the DOS has minima at around 30° and 
60° and a maximum at 45° just due to the buckling of the 
Fermi surface. To obtain the corresponding anisotropies 
of the gap, we performed strong coupling calculations for 
an Einstein phonon spectrum assuming perfect MD. We 
choose a coupling factor that reproduces in the T — M di- 
rection the experimentally reported dip structure (A w 3) 
jl4j . The absolute value of the gap depends on the con- 
sidered phonon frequency, that is why we show in figure 
4b only the relative variations of the gap. Considering 
however f2 « AOmeV as the study of the gap ratio spec- 
tral dependence [ p"7| and Raman experiments [ jTcj ] indi- 
cate, we reproduce the gap and T c of these materials 

The local minima of the gap at « 30° and 60° § re- 
flect the local minima of the DOS. Only in the case of an 
isotropic s-wave interaction in the MD regime, fine struc- 
tures of the DOS anisotropy like the local maximum at 
45°, are reflected in the gap anisotropy. If the interaction 
were not isotropic (as for example in the case of d- waves) , 
the anisotropies of the interaction would completely dom- 
inate the fine structures of the DOS anisotropies and the 
local maximum at 45° should be absent. As for the mag- 
nitude of the anisotropy, this is very dependent on the 
distance of the Van Hove singularity to the Fermi level 
. In fact the experimental uncertainty on the value of 
SE H is such that one can easily accommodate the gap 
anisotropies reported by both Refs. jl| and B, and the 
momentum dependence of the dip visibility |l4| . 

Finally, with the concept of MD we can also under- 
stand the qualitative similarity of the tunnel data and 
the ARPES data in the direction optimal for supercon- 
ductivity. In the case of MD, the tunnel spectrum is 
a sum of independent contributions from various parts 
of the Fermi surface, and it does not reflect an aver- 



aged superconducting behavior [|14j. The tunnel spectra 
are dominated by the contribution of the optimal part 
around T — M, since the Van- Hove singularity below Ep 
is extended and covers about 30% of the Brillouin zone 
and the coupling is much stronger in this region. With 
this picture we can naturally understand the asymmetry 
of the tunnel spectra of Ref. 0. In fact, the dip struc- 
ture is seen only at negative sample bias, because the 
Van Hove singularity at T — M is below the Fermi level. 
Measuring at positive sample bias, the dynamic behavior 
reflects the density of states above the Fermi level (as 
in inverse Photoemission). The presence of the dip at 
negative sample bias and its absence at positive sample 
bias Jt| , indicates that the density of states at an energy 
of the order of A above the Fermi level, is at least 30% 
smaller to that at an energy A below the Fermi level, 
and this can be easily obtained given the presence of the 
Van Hove singularity in the T — M direction. Because 
the DOS is smaller above the Fermi surface, the coupling 
is smaller and the dip is no more visible |Kj]i4|| . 
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Figure 1: Evolution of the anisotropy ratio as a func- 
tion of the characteristic range of the exchanged mo 
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menta q c . For q c < ir/\/2 it increases sharply indicating 
the onset of MD. 

Figure 2: Temperature dependence of the gap A in 
the (ir, 0) (circles) and in the (ir, ir) (triangles) directions 
for three characteristic ranges of exchanged momenta q c . 

Figure 3: Temperature dependence of the anisotropy 
ratio for q c — 7r/4 (circles), q c = tt/12 (triangles) and 
q c = 7r/20 (squares). The increase with temperature of 
this ratio is a clear indication of MD. 

Figure 4: (a): N(Ep) as a function of the angle 4> 
measured from the X — M direction for a simple second 
nearest neighbors tight binding model with a Van Hove 
singularity at lOmeV (full line), 40meV (dashed line) and 
90meV (dotted line) below the Fermi level, (b): The cor- 
responding gap anisotropy when the Van Hove singular- 
ity is at lOmeV (triangles), 40meV (circles) and 90meV 
(squares) below the Fermi level, obtained by strong cou- 
pling calculations assuming perfect MD. 
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